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Abstract: We investigate the general properties of Unified Dark Matter (UDM) scalar 
field models with Lagrangians with a non-canonical kinetic term, looking specifically for 
models that can produce a fast transition between an early Einstein-de Sitter CDM-like 
era and a later Dark Energy like phase, similarly to the barotropic fluid UDM models 
in JCAP1001(2010)014. However, while the background evolution can be very similar 
in the two cases, the perturbations are naturally adiabatic in fluid models, while in the 
scalar field case they are necessarily non-adiabatic. The new approach to building UDM 
Lagrangians proposed here allows to escape the common problem of the fine-tuning of the 
parameters which plague many UDM models. We analyse the properties of perturbations 
in our model, focusing on the the evolution of the effective speed of sound and that of the 
Jeans length. With this insight, we can set theoretical constraints on the parameters of the 
model, predicting sufficient conditions for the model to be viable. An interesting feature 
of our models is that what can be interpreted as i^de can be < — 1 without violating the 
null energy conditions. 
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1. Introduction 

In the last decade the ACDM model ||, § has emerged as the "concordance" model ||, |j] of 
our Universe. It assumes General Relativity (GR) as the correct theory of gravity, and two 
unknown components dominating the late times dynamics: i ) Cold Dark Matter (CDM), 
responsible for structure formation, ii ) a cosmological constant A making up the balance 
for a spatially flat Universe and driving the cosmic acceleration []|, ^, [?|, ^ |j| [l(J . 

Issues (still open today) with A, viz. the so-called "cosmological constant problem" 
1 11, |i~2[, prompted cosmologists to investigate alternatives to the cosmological constant, 



mainly in the form of a dynamic component dubbed Dark Energy (DE) Jl3| , 14, 16 



17, 18 1 . By now, many independent observations support both the existence of a CDM 



component and that of a separate DE j|, ||, |T^, 20, 21, 22, ||, ^, ^3], |26| . 

However, it should be recognised that, while some form of CDM is independently 
expected to exist within any modification of the Standard Model of high energy physics, 
the really compelling reason to postulate DE has been the acceleration in the cosmic 
expansion. 
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It is mainly for this reason that it is worth investigating the hypothesis that CDM 
and DE are two aspects of a single Unified Dark Matter (UDM) component, also referred 
to as "Quartessence" , see e.g. p7], |28|, p9|, I 



31, 32 



37, 38, B9L 40, 41 



42, 43, 44, 45, [16|, [17|] (a more complete list of UDM models can be found in a recent 
review 



IS]). In UDM models a single matter component must drive both the accelerated 
expansion of the Universe at late times and the formation of structures. This poses an 
interesting challenge to model building with this single component. Indeed, the accelerated 
expansion in a feature of the homogeneous and isotropic background, while the description 
of the structure formation process requires to consider inhomogeneities, at least at the first 
perturbative order. This has to be contrasted with CDM+DE models, where the CDM is 
perturbed and drives structure formation, while in general DE cannot cluster but drives 
the acceleration of the background. 

A large variety of UDM models have been investigated in the literature, mainly based 
on adiabatic fluids or on scalar field Lagrangians. For example, the generalised Chaplygin 
gas [27, 28, 29, 31 1, the Scherrer and generalised Scherrer solutions p4|] , the perfect fluid 
with "affine" equation of state [37, |39| (see |38j for the corresponding scalar field model), 
and the homogeneous scalar field deduced from the galactic halo space- time |4|, |§. 

It is also possible to reinterpret UDM models based on a scalar field Lagrangian with 
a non-canonical kinetic term in terms of non- adiabatic fluids fH |l| (see also @, ||, 



)• 



For these models, we can use Lagrangians with a non-canonical kinetic term (namely a 
term which is an arbitrary function of the square of the time derivative of the scalar field, 
in the homogeneous and isotropic background). 

Originally, these Lagrangians were proposed to obtain inflationary solutions driven 
by kinetic energy, the so-called /c-inflation (5^, [53|. This scenario was also adopted as a 
description of DE [54, [55], 56]. When /c-inflation was extended to embrace more general 
Lagrangians [pTl |5q| , the resulting new scenario was then dubbed ^-essence (see also p9|, 
6C, 56, 62, |63|). An important feature of scalar field models is the effective speed of 
sound, which remains defined in the context of linear perturbation theory, is not the same 
as the adiabatic speed of sound (see [p4l 65, [53], 



cf. 



Independently of the approach chosen to build a UDM model with a single matter 
component, in general the latter will have pressure perturbation in the rest frame (which is 
then gauge invariant), which implies a non vanishing effective speed of sound (see j67j, cf. 



also [68, 69, 64, [6^]). The latter corresponds to a Jeans length (i.e. a sound horizon) below 
which clustering is suppressed [j67| , |35| , |39| . Moreover, the evolution of the gravitational 
potential on scales smaller than the Jeans length is characterised by oscillations and decay, 
thus producing a strong late-times Integrated Sachs Wolfe (ISW) effect which may spoil 
the agreement with the observed Cosmic Microwave Background (CMB) radiation angular 
power spectrum [35]. Therefore, the general lesson to be drawn is that a viable UDM model 
must necessarily be characterised by a vanishingly small effective speed of sound, so that 
structure formation would not be suppressed (at least in the linear regime) and the ISW 
effect would be compatible with CMB observation. 

However, for many UDM models the requirement of a vanishingly small speed of sound 
implies such a severe fine-tuning on the parameters that they become practically indistin- 
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guishable from ACDM, see for example [31, 32, 33, 39, 7C]. Fortunately, this does not need 



to be the case in general: indeed, in order to avoid such undesirable feature, some possible 
way-outs have recently been proposed. In we introduced a class of UDM models whose 
equation of state allows for a fast transition between an early CDM-like era and a later 
ACDM-like phase. The speed of sound can be very large during the transition, but if the 
latter is fast enough, these models predict a large-scale structure spectrum and an angular 
CMB temperature spectrum in agreement with observation. Ref. |^5|] explored unification 
of DM and DE in a theory containing a scalar field of non-Lagrangian type, obtained by 
direct insertion of a kinetic term into the energy-momentum tensor. Ref. |^7j introduced a 
class of models based on two scalar fields, one of which is a Lagrange multiplier enforcing 
a constraint between the other scalar field and its derivative. The purpose is to assure 
a vanishing speed of sound. In [Q, the authors devised a reconstruction technique for 
Lagrangian which allows to find models where the effective speed of sound is so small that 
the ^-essence scalar field can cluster (see also pi, ||, gf|). 



In the present paper we develop the technique of [|4(j] to build scalar field Lagrangians 
with non-canonical kinetic term which can produce a fast transition, similarly to the 
barotropic fluid UDM models we proposed in [44|. However, while the background evo- 



lution can be very similar in the two cases, the perturbations are naturally adiabatic in 
fluid models, while in the scalar field case they are necessarily non-adiabatic |34], ^5|, |66| , 
cf. n, H HI- This 

new approach allows us to escape the problem of the fine-tuning on 



the parameters which usually plague many UDM models. 

The rest of the paper is organised as follows. In Section || we present the basic equa- 
tions describing the background and the perturbative evolution of a general /c-essence 
Lagrangian. In Section ||| we generalise the class of UDM models investigated in |^] and in 
Section ||] we introduce a UDM model with fast transition, analysing its background evolu- 
tion. In Section [| we analyse the properties of perturbations in this model, focusing on the 
the evolution of the effective speed of sound, of the Jeans length and of the gravitational 
potential during the transition. Section ^ is devoted to our conclusions. 
Throughout the paper we use 8ttG = c 2 = 1 units and the (—,+,+,+) signature for the 
metric. Greek indices run over {0, 1, 2, 3}, denoting space-time coordinates, whereas Latin 
indices run over {1,2,3}, labelling spatial coordinates. 

2. Background and perturbative equations for UDM scalar field models 

The action describing a generic Unified Dark Matter scalar field model in GR can be written 
as 

' ' - (2.1) 



S = Sg + S v = j d x^J^g 

where 



1 

2 

The energy-momentum tensor of the scalar field (p is 



X = --V^VV • (2-2) 



= -77^f£ = ^ (iy «V„ y - + IX,. Xu hll , (2.3. 
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If Vfj,(f is time-like, S v describes a perfect fluid with Tfi v = (p + p)u^u u + pg^ u , with 
pressure 

C = p(tp,X), (2.4) 



and energy density 



p = p(<p, X) = 2X 9p ^ x X) - p(<p, X) , (2.5) 



where 

(2.6) 



r V2X 

is the "energy frame" four- velocity [71, In this frame, the kinetic term becomes X 



(p /2, where (p = u^V^ is the proper time derivative along u^. 

Now we assume a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) background 
metric with scale factor a(t). When the energy density of radiation becomes negligible, and 
disregarding also the small baryonic component, the background evolution of the Universe 
is completely described by the following equations: 

r2 = \pi ( 2 - 7 ) 
H = ~(p + p), (2.8) 

where the dot denotes differentiation with respect to the cosmic time t (coinciding with 
proper time along u^) and H = a/a is the Hubble parameter. Eqs. fl2.7j )- (|2.8| ) imply the 
energy conservation equation 

P = -3H(p + p). (2.9) 

On the background FLRW metric, the equation of motion for the scalar field (p(t) 
follows from the energy momentum tensor (2.3) and Eq. 



The two relevant quantities that characterize the background and perturbative dynamics 
of a UDM model are the equation of state w = p/p and the effective speed of sound c s , 
which, in our case, read 

p 

w = 2X(dp/dx)- P ' {2 - n) 
2 _ dp/dx _ dp/dx 

s ~ dp/dx (dp/dx) + 2X(d 2 P /dx 2 ) ' 1 ' ' 

The effective speed of sound plays a major role in the evolution of the scalar field per- 
turbations 5(p and in the growth of the overdensities 5p. In fact, let us consider small 
inhomogeneities of the scalar field, cp(t,x) = <po(t) + 5<p(t,x), and write FLRW metric in 
the longitudinal gauge: 

ds 2 = -(1 + 2<$>)dt 2 + a 2 (t)(l - 2<$>)5 ij dx i dx j , (2.13) 
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where Sij is the Kronecker symbol and we have used the fact that 8T- = for i ^ j [ff2fl. 



Indeed, for this perfect fluid case, there is a unique peculiar gravitational potential <I>, 
analogous to the Newtonian potential [71]. When we linearize the (0 — 0) and (0 — i) 
components of Einstein equations (see Refs. [53, |66|]), we obtain the following second order 
differential equation [66, 35]: 

_ _ ft" 

(2.14) 



c 2 V 2 u - —u 







where u = 2$ / {p + p) 1 / 2 , 6 = (1 + p/ p)~ l l 2 /{y/la) and the prime denotes derivation w.r.t. 
the conformal time r/, which is defined through drj = dt/a. As it is clear from ( p.lOj ) and 
J2.14 ) , the effective speed of sound plays a crucial role in the evolution of the scalar field 
perturbations Sip and of the gravitational potential. 



3. Prescriptions for building a scalar field UDM model 

A well-known issue of UDM models is the existence of an effective sound speed [see 
Eq. (|IJ)], which may become significantly different from zero during the Universe evo- 



lution and thus leading to the appearance of a Jeans length (i.e. a sound horizon) below 
which the dark fluid cannot cluster (e.g. see |37], ^]). Indeed, the presence of a non- 
vanishing speed of sound modifies the evolution of the gravitational potential, producing 
a strong Integrated Sachs Wolfe (ISW) effect [35]. For this reason, one of the main issues 



in the framework of UDM model building is to investigate whether the single dark fluid is 
able to cluster and produce the cosmic structures we observe. 

In Ref. [4C] the authors proposed a technique to construct UDM models where the 



scalar field mimics the ACDM background evolution and, at the same time, has an effective 
sound speed that is small enough to allow for structure formation, also avoiding a strong 
ISW effect (see also ©EH). 

In this Section we develop and generalize the approach of Ref. p0|] . Specifically, we 
focus on scalar field Lagrangians with non-canonical kinetic term to obtain Unified Dark 
Matter models where a single component can mimic the dynamical effects of Dark Matter 
and Dark Energy. This sets the ground for Section [|, where our aim is to build scalar 
field models in which the dynamics is characterised by a fast transition between an early 
CDM-like phase, with the background following an Einstein de Sitter evolution, and a 
late accelerated DE-like phase, in this way generalising the UDM fluid models with fast 



transitions we presented in [70]. 



We look for a scalar field Lagrangian £ whose classical trajectories are directly de- 
scribed by an appropriate pressure p(N) that we choose a priori, where N = In a. Let us 
consider a Lagrangian of the form 

£(<p, X) = p(ip, X) = f{ip)g{h{ip)X) - V{ip) . (3.1) 

Note that, with the freedom allowed by the three potentials f(<p), h(ip) and V(tp), we can 
independently specify the equation of state parameter w and the sound speed c s ^0|. In 
order to reconstruct these potentials we need three dynamical conditions: a) a choice for 
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p(N), b) the continuity equation or, equivalently, the equation of motion ( 2,10 ), c) a choice 
for c 2 s (N). 

Now we look at the four main prescriptions needed to to obtain the various terms of 
the Lagrangian (3.1). 



3.1 Equation of state of UDM 

Assuming p(N) as a function that we choose a priori, Eq. 

dp(N) 



becomes 



i.e. 



p(N) = e 



dN 



-3N 



+ 3p(N) = -3p(N) , 



/N 
(e 3N 'p(N>)dN') + K 



(3.2) 



(3.3) 



where K is an integration constant. In particular, by imposing the condition 
C(X, ip) = p(N) along the classical trajectories on cosmological scales, we obtain 



M 



p(N) 



Here, we are constraining the solutions of the equation of 



motion to live on a particular manifold A4 p rm [depending on our choice of p(N)] embedded 
in the four dimensional space-time. 

We assume that the pressure p and energy density p of our UDM: a) satisfy the null 
energy condition p + p > [73, [5^] (here we are not interested in phantom cosmology); b) 
violate the strong condition at late time, i.e. p < —p/3, in order to reproduce the observed 



accelerated expansion of the Universe; c) imply an attracting fixed point for Eq. fl3.2|) : for 
N — > oo, p — > pa and p — >■ —pa, so that dp/dN — > and pa plays the role of an unavoidable 
effective cosmological constant [j74|, |75|, [37| . 

Now, if p = — pa, then K = (po — Pa), where po is the value of p today. The freedom 
provided by the choice of K is particularly relevant. In fact, by setting K = 0, we can 
remove the term p oc a -3 from Eq. ( |3.3| ). Alternatively, when K ^ 0, we always have a 
term that behaves like pressure- less matter (see Ref. pQ|), with density today K = p m Q = 
(po — Pa)- This simply follows from assuming p(N) as given, so that, in the solution p(N) 
J3.3D , Pm{N) = p m (>e~ 3N represents the matter-like homogeneous solution of (|3.2|) , and 



p(N) generates the particular solution (see [48|, [40], |7q|). 

We conclude with the following two remarks: i) starting from Eqs. (3.2) and ( ft.3j ) we 
can obtain a dark component which can mimic not only a cosmological constant but also 
any dynamical DE; ii) in order to obtain explicitly the three potentials f(f), h(ip) and 
V((p), p(N) should be chosen in order to have an analytic expression for p(N). 



3.2 Reconstruction of the potential f{(p) 

From the Lagrangian C = f((p)g(h((p)X) — V(tp), and Eq. (p^), we get 
■dg(h(<p(X,N))X)] tr _ rv „ m __, An , _ 3J v r fV 



2X 



dX 



f(cp(X,N))=p(N) + e 



/I V 
(e 37V V(A r ')^ / ) + K 



(3.4) 
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The relations (|3^ ) and tp = & 1 \X,p(N))\ M enable us to determine a connection 



p(N) 



between the scale factor a and the kinetic term X on the manifold M p [n) and, therefore, 
a direct mapping between the cosmic time and the manifold M. p ij^\- 
Let us define f(<p(X,N)) in the following way 



f(v(X,N)) 



p(N) + p{N) 



2X [dg(h{<p{X,N))X)/dX] 



A(N,X) , 



(3.5) 



along the classical trajectories on cosmological scales, we deduce that A(N, X) = 1. Specifi- 
cally, let us define A(N, X) = A(N, X/p(N)). Then, knowing that X = (dN / dt) 2 (dp / dN) 2 /2 
p{N)(dp/dN) 2 /6, we get 



X = p(N)A(-V{N) , 

f N / ~ \ 1/2 

= v»i± / (sA^iN')) <W . 



(3.6) 
(3.7) 



Without any loss of generality, consider the case with the + sign in front of the integral 
above. Redefining Eqs. ( [3.6D and ( [3.7D as X = Q P {N) and p = Q P (N), we can write 
f(X,N) = f(g p (N),N) = f@ P {<^b)),& X \<p)) = f(<p). Obviously, in order to 
obtain f(p) we must have chosen an appropriate g(h(<p)X) (see point 3). 
Now, let us assume A(N,X) being of the following form: 



A(N(a),X) 



Pa 3 (ip(N(a))/(2Pv) -X 



X 



(3.8) 



where (3 and p, are two appropriate constants whose choice depends on p(N) and c 2 (N) and 
v = S7 m /r2DE 5 where ri m = p m o/(3ifg) and Ode = 1 — ^m- Along the classical trajectories 
on cosmological scales, we obtain 



X(a) = G p (N(a)) : 
p(a) = Q P (N(a)) 



V p(N(a)) 
1 + W/3)a~ 3 

6vp 



— ) a2 



(3.9) 
(3.10) 



In conclusion, from Eq. (3.5), relations ( |3.9| ) and (3.10) allow us to determine explicitly 
f(ip). Obviously if p = —p\ we get the same X(a) and ip(a) obtained in Ref. p0|]. 

3.3 Reconstruction of the potentials h(p) and V(<p) 

In order to get h(p) we have to know g(X), where X = h(p)X. From ( [2. 12 ) we have 

dg(X)/dX 



(dg(X)/dX) + 2X(d 2 g{X)/dX 2 ) 



(3.11) 



where we must impose (? S {X) > 0. In general c 2 is not bounded from above; however in 
specific cases the condition c 2 (X) < 1 must be satisfied. Knowing that X = X(a) from 
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Eqs. Q and ( g7lO|) , we can choose an appropriate function c 2 (X) = 1Z p (N(a)) such that 
(? s -C 1. Therefore, 

= — (3.i 2 ) 

H p{N{a)) 

and, through Q p ^ (</?), we can reconstruct /i(</?)- 
In order to get V(<p) we have to use Eq. fl3,l|): 



V(p(a)) = /(^(a)) 5 (^(a)) - p(N(a)) . (3.13) 
Using Qp ((f), we can reconstruct F((/?)- 
3.4 Simplification of the Lagrangian 

At this point, it is important to stress that the Lagrangian we have obtained may be very 
complicated. In [^] it was shown that there exists a class of Lagrangians having similar 
kinematical properties: in particular the equation of state w and speed of sound c s are 
invariant under certain transformations within the class. Thanks to this property we are 
able to remarkably simplify our Lagrangian. 

Specifically, from C(X, <p), we can obtain a new Lagrangian C(R(cp)Y, (p) where 



2 R(cf>(<p)) 



and <t>(<p) = ± [R(fiy 1/2 dp] + K , (3.14) 



with R((f>) > and where K is an appropriate integration constant. Without any loss of 
generality, consider the case with the + sign in front of the above integral. Defining 



[R{4>((p))]~ l l 2 = cosh 
we obtain the following relations 

4/3 



va 3 



3^/3 V' 2 



1/2 



(3.15) 



(3.16) 



R{m) = i + pw(w ■ (3 - 17) 

3.5 Prescriptions: further remarks 

Let us stress that the receipt described in the previous subsections is completely general 
for any p(N) (when p + p > 0) and for any positive c 2 (N). Obviously, the values of c 2 (N) 
depend strongly on the definition of g(X), through Eq. ( f3.ll ). 
Now additional comments are in order: 

i) From Section [3.1| we note that, having assumed p m (a) = p m Q a 3 as the pressure-less 
solution of Q3.2j ), we can formally define 

wde = — - — , (3.18) 

P- Pra 
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where p — p m is the DE-like part on our UDM energy density and p is its pressure. 
Interestingly, it is possible to build models with u>de < — 1 without violating the null 
energy condition, see the next Section for explicit example. 

ii) Choosing an arbitrary equation of state w(N) we obtain 

p(N) = pi w ( N)e -3f N (HN')+l)dN> ; (3 _ 20) 



where pi is a positive integration constant. Therefore, imposing again the condition 
C(X, (p) = p[w(N), N] along the classical trajectories, i.e. ip = 
we get 

2X MK<p[x,m)x) f{XjN) = pi [w{N) + l]e -3f N ( W (N') + i) d N* (3 21) 

oX 

Hence, on the classical trajectory we can impose, via w(N), a suitable function p(N) 
and thus the function p(N). Also in this case, from Eq. (|3 .21 ) and following ar- 
guments similar to those described in the points above, we can get the relations 
X = Q W (N), and consequently 



M w(N y 



V=Q W (N) = ± f l[6G w (N')} 1/2 



p . e -3f N ( W (N")+l)dN" 



-1/2 ^| 

dN' \+ipi. 



(3.22) 

With the functions Q W {N) and Q W (N), we can write f(X,N) = f(g w (N),N) = 
f{Gw(Qw 1 ^(v 9 ))i Qw = /(^)- Then we can find a Lagrangian whose behaviour 

is determined by w(N) and whose speed of sound is determined by the appropriate 
choice of g{h{ip)X). 

iii) Once we have obtained a Lagrangian, it is important to investigate the kinematic 
behaviour of the UDM fluid during the radiation-dominated epoch and, in particular, 
to choose appropriate initial conditions in order to ensure that, during recombination, 
we obtain solutions of the equation of motion which properly describe p(N) during 
the subsequent dark epoch. 

In the next section, we give an example in which we apply the prescriptions above and 
consider the explicit solutions in the case of the following Dirac-Born-Infeld type kinetic 
term: 



g{h{ip)X) = - v / l-2h(^)X , (3.23) 

which implies 

c 2 s (h((p)X) = 1 - 2h((p)X . (3.24) 

Moreover, we will choose a suitable p(N) in order to obtain a class of UDM models charac- 
terized by a fast transition in the equation of state and, in section [|, an appropriate c^(N) 
such that the /c-essence scalar field can cluster. 
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4. A simple UDM scalar field model with fast transition. 



We now propose a UDM model with the following equation of state, which we give in 
parametric form: 



p(N(a)) 



PA 
2 



1 + tanh 



(a 3 - a 3 ) 



f AX(u)) - ,, A { ~ + A a 3 In ^ cosh 



(a 3 - a 3 ) 



+ PmO a - 3 
PA 



(4.1) 
(4.2) 



where (|4.2| ) follows from (44.) by integrating Eq. (|3.3| ). The parameters ot and (3 respec- 
tively represent the scale factor value at which the transition takes place and the rapidity 
of the transition. The third parameter of the model is p\: it is clear from Eq. ( |4.2|) that in 
the limit a — > oo p — > p\ and p — > —p\, so that p\ plays the role of an effective cosmological 
constant, i.e. a fixed point of Eq. ( |3.2[ ) (see Section |3.1[ ). Finally, we note that our UDM 
model depends on a fourth parameter, p m o/pA = ^m/^A in Eq. ( p~2j ) = pa/{3Hq)). 
In comparison with our barotropic fluid model in JTQ| , here we have this extra parameter 
because we have explicitly introduced a matter-like part in the energy density thanks to 
the integration constant K in (|3.3| ). 
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Figure 1: Illustrative plot of p/pa, p/pa and w as functions of the redshift z, with zt = 2 and 
^m/^A = 3/7. The lines, from short to long dashes, correspond to (3 = 1, 10, 100, respectively; the 
black solid line corresponds to f3 = 1000. For reference, in the bottom panel, we also plot the total 
w for the ACDM model (thin black line) with Sl m /Q A = 3/7 (cf. @ ^ |i[). 
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In Fig. |H we show the evolution of p(N(z)), p(N(a)) and the equation state parameter 
w(N(a)) as a function of the redshift z. An important difference with the barotropic 
model in JFTJ is that the parametric representation (4.1)- (4. 2) is effectively equivalent to 
an equation of state p = p(p,s), where s is an entropy density. Therefore, in general the 
condition p = —p does not uniquely determine a fixed point of Eq. (3.3) as in the barotropic 
case. In other words, in a certain region of parameter space the condition p = —p has two 
solutions: one is the effective cosmological constant p\, the asymptotic value of p and —p 
for a — > oo, the other is a minimum value of energy density < p\ that is attained in 
the future for a finite value of a* > 1, with p{a*) = — /?*. For a > a* the null energy 
condition is violated, so that the equation of state becomes phantom and p grows again, 
asymptotically approaching p\. In this paper we are interested in mapping the equation 
of state (4.1)-(fi~2]) into a Lagrangian for a scalar field that does not violate the null energy 
condition, thus we will focus on the relevant region in parameter space (see below). Fig. |2|is 




5 10 15 20 25 

P 

Pa 



Figure 2: Illustrative parametric plot oip/pA as a function of p/pa, with Zt = 2 and il m /rt\ = 3/7. 
The short to long dashed lines correspond to (3 — 1, 10, 100, respectively; the black solid line 
corresponds to (3 — 1000. For reference we also plot the p = —p line. Note that the latter line 
is crossed by the f3 = 1 model (see text). All models asymptotically evolve toward the effective 
cosmological constant p\. 



a parametric plot of p{a) vs. p(a), assuming a transition at zt = 2 and for a representative 
choice of the other parameters. The curve for (3 = 1 illustrate in particular a case where 
indeed the parametric equation of state ( f4.1|) -(|4~2^) becomes phantom in the future. 

The parametrization of Eqs. (|4.1| )- (|4.2| ) in terms of p\ is mathematically natural, but 
it can be deduced from Figs. ||-|| that it is not so practical from a phenomenological point 
of view; a more useful parametrization is obtained using pde = Po — PmO, he. the present 
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value of the DE-like part of our UDM. We obtain the following relation: 



A ^ + hi <i cosh 



§(i-8 



(4.3) 



Pde 2 2/3 

In Fig. || we plot pde/pa for different values of (3 and z t . Notice that, for z t > 2 and large 



i. 
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Figure 3: Left panel: pde/pa as function of z t , for j3 — 0.1, 1, 10, 10 2 , 10 3 (from bottom to top). 
Right panel: pde/pa as function of j3, for z t — 0.5, 1, 1.5, 2, 2.5, 5 (from bottom to top). 



(3, we have that pde/pa — > 1- 
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Figure 4: The gray region represents the part of the Zt-/3 plane where the null energy condition 
p + p > is satisfied, assuming the representative value v = 3/7. 
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Figure 5: Illustrative plot of p/puE, p/pde and w as functions of the redshift z, with zt = 2 and 
v = 3/7. The lines, from short to long dashes, correspond to /3 = 4, 10, 100, respectively; the black 
solid line corresponds to /3 = 1000. For reference, in the bottom panel, we also plot the total w for 
the ACDM model (thin black line) with fi m /0 A = 3/7 (cf. @ ^ |§). 



Now, combining Eq. ( ^4.3[) with Eqs. (4.1)- fl4,2|) , p(a) and p(a) can be re-expressed as 



(l + tanh[(/3/3) (a 3 -a t 3 )]} 
" ^ 1 + 0//3) In {cosh [03/3) (l - a 3 )] } ' ^ 

f l + (3//3)a- 3 ln{cosh[(/3/3)(a 3 -a t 3 )]} _ 3 \ 
= { l^^lnlcosh^a-a 3 )]} + ^ ) ' 



where v = p m o/pDE = fi m /^DE (see Section |3.2| ). A good guess for the value of this 
parameter is v = 3/7: assuming this value, in Fig. ||| we plot the region of the Zt — /3 plane 
where the null energy condition p + p > is satisfied. 

In Fig. we plot p(z), p{z), w(z), as functions of the redshift z for different values of 
j3 and z t and for v = 3/7. 

Finally, in Fig. ^ we show a parametric plot of p(a) vs. p(a), assuming a transition at 
z t = 2 and for v = 3/7. From Figs. [7| and |8] and Eq. ( |3.18| ) we see that for z < z t we can 
obtain u>de < — 1 without violating the null energy condition [see comment i) in Section 



3.5 1- When (3 > 10 3 and z t > 2, we have that ^de — PA] in this case, our class of UDM 
models is characterised by a fast transition regime and today u>de — — 1- However, this is 
not the case for (3 < 10 3 or z t < 2, so that wbe can be significantly smaller than —1 even 
at small redshift. 
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Figure 6: Illustrative parametric plot of p/pde as a function of p/puE, with zt = 2 and f = 3/7. 
The short to long dashed lines correspond to (3 — 4, 10, 100, respectively; the black solid line 
corresponds to j3 = 1000. 
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Figure 7: Illustrative plots of (p + p)(z)//°de an d wde as a function of the redshift z, with 
zt = 2 and v = 3/7. In both panels the short to long dashed lines correspond to j3 = 4, 10, 100, 
respectively; the black solid line corresponds to = 1000. Clearly (p + p) > always, i.e. the null 
energy condition is never violated. 



5. Analysis of the Jeans wave number and the gravitational potential 

Knowing p(N) and p(N), following the prescriptions of section || we still need to choose a 
suitable (? a in order to obtain the Lagrangian C(Y, <ft) that will completely specify our UDM 
scalar field model. This Lagrangian will reproduce our choice of p(N), p(N) and c%(N) 
on-shell, i.e. along the classical trajectories on cosmological scales (see Appendix [A]). In 



this section we choose a suitable cj and, via the equation state (|4.4j) -( [4.5D , we study the 



evolution of the Jeans wave-number and the gravitational potential. 
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Figure 8: u>de as a function of the redshift z, with v = 3/7 and z t — 1 (left panel) and z t = 6 (right 
panel). In both panels the short to long dashed lines correspond to j3 = 5, 10, 100, respectively; the 
black solid line corresponds to j3 = 1000. Without violating the null energy condition our models 
produce wde after the transition. 



5.1 Speed of sound and Jeans scale 

We assume a speed of sound of the following form: 

„2 



1 + (1/0-3/4,)" ' 



(5-1) 



where Coo > and n > are free parameters, the former representing the asymptotic future 
speed of sound. 

In Figs. H] and |l(] we have plotted cf for different values of Cqo and n. We can deduce 
that, for practicality and in order to reduce the number of free parameters, we can choose 
n = 4, 5 and Cqo = 0.1. In this case we obtain a value for the speed of sound that is suitably 
small at all times, without requiring any fine tuning. 

Let us now analyse the Jeans scale; rather then directly considering the physical Jeans 
length A j we look at the Jeans wave number kj = 2Tra/\j. Starting from Eq. ( 2.14 ), the 
squared Jeans wave-number is defined as follows [j35f: 



(5.2) 



This is a crucial quantity in determining the viability of a UDM model, because of its 
effect on perturbations, which is then revealed in observables such as the CMB and matter 
power spectra. Indeed, any UDM model should satisfy the condition kj 3> k 2 for all the 
scales of cosmological interest, in turn giving an evolution for the gravitational potential 
<&(?7, x) oc $(?7, k) exp (ik • x) of the following type: 



a(rj) 



a(rj) 2 dfj 



(5.3) 



where A k = (0, k) T m (k), $ (0, k) is the primordial gravitational potential at large scales, 



set during inflation, and T m (k) is the present time matter transfer function, see e.g. [77] 
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Figure 9: Illustrative plot of c 2 as a function of the redsliift z, with v = 3/7 and n = 2 (left 
panel) and n — 4 (right panel). Left panel: the short to long dashed lines respectively correspond 
to Coo = 10 _1 , 10~ 2 , 10~ 3 and the black solid line corresponds to Coc = 10~ 4 . Right panel: the 
short to long dashed lines respectively correspond to Coc = 0.5,10 _1 ,10 -2 and the black solid line 
corresponds to Coc = 10 ~ 3 . 
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Figure 10: Illustrative plot of c 2 as a function of the redshift z, with Coo = 0.1 and v — 3/7. 
The lines, from short to long dashes, respectively correspond to n = 2, 4, 6 and the black solid line 
corresponds to n = 8. 



From (|T^) , the explicit form of the Jeans wave number for a scalar field UDM model 
turns out to be 



A: 



P 



(l + w) 



2 (1 + z) 2 c 2 s 



2l C ad 



W 



4)' 



+ 



6(1 + w) 



1 

+ 3 



(5.4) 



where c 2 d = p'/p'. Comparing with the Jeans wave number for a barotropic fluid (see Eq. 
(3.3) in [44|), for which c 2 = c 2 d , the only difference is the overall 1/c 2 factor replacing 
l/c 2 d : this gives extra freedom in building a suitable kj. 

Clearly, we can indeed obtain a large /c 2 when c 2 — > 0; in addition, when c 2 d changes 
rapidly around zt, i.e. when the above expression is dominated by the \p (c 2 ^)'// 9 '] term, 
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we can have a fast trasition. Therefore, we can conclude that, for (3 < 1000 and before and 
after the fast transition for (3 > 1000, it is crucial that Cg be sufficiently small. Defining 



Cg as in Eq. ( |5.1|) , in Figs. 11 we plot the Jeans wave- number for the representative case 



Zf 



2 and Cqo = 0.1, for various values of /3 and n. 



£,(h-Mpc ' 



k,(h-Mpc~ 




Figure 11: The Jeans wave-number fcj (in h Mpc 1 units) assuming v — 3/7, z% — 2, Coo = 0.1 
and P = 4 (left-top panel), /3 = 12.1 (right-top panel), (3 = 100 (left-bottom panel) and (3 = 1000 
(right-bottom panel). The lines, from short to long dashes, respectively correspond to n = 3,4, 5 
and the black solid line corresponds to n = 6. 



From the right panels of Fig. 11, we can note that for j3 = 12.1 and 1000 the Jeans 
wave number kj momentarily vanishes. In general, around these points the corresponding 
Jeans length becomes very large, possibly causing all sort of problems to perturbations, 
with effects on structure formation in the UDM model. On the other hand, for sufficiently 
small c s we note that i) in general the Jeans wave number becomes larger and ii) it becomes 
vanishingly small for extremely short times, so that the effects caused by its vanishing are 
sufficiently negligible, as we are going to show in the next subsection when we will analyze 
the gravitational potential <F Therefore, in building a phenomenological model, we can 
choose its parameter values in order to always satisfy the condition k <C fcj for all k of 
cosmological interests to which linear theory applies. In other words, we can always build 
our model in such a way that all scales smaller than the Jeans length A -C Aj correspond to 
those in the non-linear regime, i.e. scales beyond the range of applicability of the model. So, 
for these scales, no conclusions on the behaviour of perturbations can be derived from the 
linear theory. Indeed, to investigate these scales, one needs to go beyond the perturbative 
regime investigated here, possibly also increasing the sophistication of the UDM model in 
order to properly take into account the greater complexity of small scale non-linear physics. 
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Moreover, from Fig. 11 we can finally conclude that, for n = 4 or 5 and c c 
a acceptable Jeans wave number at all times. 



0.1, we obtain 



5.2 The gravitational potential 

From Eq. Q2.14 ) let us write explicitly the differential equation for the gravitational poten- 
tial <3?: 

d 2 ${a,k) 



da 2 



i 1 dn i 4 1 z c *A d ^ k) | 

H da a a I da 



2 dTi 1 M 2 , c 2 k 2 
aH da a 2 1 



a 2 n 2 



*(a,fc) = , 
(5.5) 

where % = a 2 H 2 and plane-wave perturbation, $(a,x) oc <l>(a, fc) exp (ik • x) have been 
assumed. 



In Figs. 12 we plot the gravitational potential for different values of n, f3 = 0.1, zt = 2, 
Coo = 0.1 and for k = 0.05 /i Mpc^ 1 and k = 0.2 h Mpc" 1 . We note that, for k ~ fcj, 
&(a,k) starts to oscillate and decays, thus preventing structure formation. On the other 
hand for n > 2 and for small values of /3, we obtain a shape of the gravitational potential 
very close to that of the ACDM model. 




Figure 12: Gravitational potential <£>(z, k) as a function of the redshift z, assuming v ~ 3/7, z% = 2, 
Coo = 0.1 and j3 — 4. The lines, from short to long dashes, respectively correspond to n = 1, 1.5,2. 
For comparison, the black solid line represents $ the ACDM model with J7 m /f2A = 3/7. Left panel: 
k = 0.05 h Mpc -1 ; right panel: k = 0.2 h Mpc^ 1 . 



Now, assuming for simplicity n = 4 and Coo = 0.1 (in order to have an acceptable fcj) 
we analyse the gravitational potential and investigate how it depends on the background 
parameters /3 and z t (or, equivalently, at). As we know from Fig. |5[ for j3 < 100 the value 
of Zt practically loses the meaning of scale parameter at the transition. For these range of 
parameters, we can also observe this effect in Fig. |l^, where we have plotted $>(z, k) for 
different values of (3 and zt- 

Moreover, from these panels we can note another interesting effect. For large values 
of /3 and for small values of zt, &(z,k) is constant in time (as it should be in a pure 
matter Einstein De Sitter model) and assumes the same value of $(0, k) until z ~ Zt 
and then for z < z t its value quickly goes down, eventually intersecting the potential of 
the ACDM model. This happens because the background starts to "feel" the effective 
cosmological constant only for z < z t . It is important to stress that this evolution of $ 
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Figure 13: Illustrative plots of the gravitational potential 4>(z, k) as a function of the redshift z, 
for k = 0.2 h Mpc -1 , v — 3/7, z t — 1, Coo = 0.1 and n = 4. For comparison, the black solid line 
represents $ in the ACDM model with il m /Cl\ = 3/7. Left-top panel: z t = 1. The lines, from 
short to long dashes, respectively correspond to j3 = 10 4 , 50, 5. Right-top panel: z t — 2. The lines, 
from short to long dashes, respectively correspond to (3 — 10 4 ,50,4. Left-bottom panel: z t = 3. 
The lines, from short to long dashes, respectively correspond to /3 = 10 4 , 50, 4. Right-bottom panel: 
Zt = 5. The lines, from short to long dashes, respectively correspond to /3 = 10 4 ,50,4. 



could produce a strong Integrated Sachs Wolfe (ISW) effect: however, this would only be 
due to the particular evolution and would not depend from d?. Obviously, this effect is 
weaker if we set zt > 2 and completely negligible for zt > 5, i.e. these UDM models become 
indistinguishable from the ACDM model, cf. 44]. 

In conclusion, the suitability of models with 4 < (3 < 1000, as well as that of models 
with (3 > 1000 and Zt < 2, need a study of the matter and CMB power spectra, which 
will reserve for the future. On the other hand, for (3 > 1000 and assuming z t > 2, i.e. for 
an early enough fast transition, the above analysis shows that our UDM models should be 
compatible with observations. 

In order to compare the predictions of our UDM model with observational data, we 
have to define the UDM density contrast as 5 := 6p/pA [p9Q , where here pa = P — PA 
is the clustering "aether" part of the UDM component 0, p6| . Indeed, our equation of 
state admits an asymptotic (a — > oo) effective "cosmological constant" [[74|] which we have 
already defined as pa in section ||. In this case, starting from the perturbation theory 
that we outlined in section ||, we can infer the link between the density contrast and the 
gravitational potential via the Poisson equation for scales smaller than the cosmological 
horizon and z < z Tec , where z Tec is the recombination redshift (z rec ~ 10 3 ) in the following 
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way: 



5(k,z) 



2k 2 $>(z,k) (1 + z) 
PA 



(5.6) 



6. Conclusions 



H|, ^3fl , the search 

[/]. |l0| and the measurements of the CMB anisotropics 
25] are very well explained by assuming that two dark components govern the 



The last decade of observations of large scale structure [20, 21, 22 
for type la supernovae (SNIa 

1- E 



dynamics of the Universe. They are DM, thought to be the main responsible for structure 
formation, and an additional DE component that is supposed to drive the measured cosmic 



acceleration [17, [T^, 16, 14, 15]. However, it should be recognised that, while some form of 



CDM is independently expected to exist within any modification of the Standard Model of 
high energy physics, the really compelling reason to postulate DE has been the acceleration 
in the cosmic expansion. It is mainly for this reason that it is worth investigating the 



hypothesis that CDM and DE are two aspects of a single UDM component, see e.g. [48, 18]. 

In this paper we have developed and generalized the technique to construct scalar field 
UDM models proposed in Ref. ]j40| and we have focused on Lagrangians with non-canonical 
kinetic term to obtain models where a single component can mimic the dynamical effects 
of Dark Matter and Dark Energy and, at the same time, has a sound speed small enough 
to allow for structure formation. 

In the second part of the paper, we have built UDM models which can produce a fast 



transition, similarly to the barotropic fluid UDM models we proposed in [44]. However, 
while the background evolution can be very similar in the two cases, the perturbations 
are naturally adiabatic in fluid models, while in the scalar field case they are necessarily 
non-adiabatic 



, 65, 66], cf. [67, 38 



|. This new scalar field model allows to escape the 
problem of the fine-tuning on the parameters which usually appeared in many previous 
UDM Lagrangians. 

First of all, an interesting feature of our models is that, for z < Zt, u>de can De < — 1 



without violating the null energy conditions [see comment i) in Section 3.5]. Subsequently, 
we have analysed the properties of perturbations in our model, focusing on the evolution 
of the effective speed of sound and that of the Jeans scale. In general, in building a 
phenomenological model, we have chosen its parameter values in order to always satisfy 
the condition k <C kj for all k of cosmological interests to which linear theory applies. In 
this way, we have been able to set theoretical constraints on the parameters of the model, 
predicting sufficient conditions for the model to be viable. In particular, we have found 
that for sufficiently small c s , i) the Jeans wave number becomes larger and ii) it becomes 
vanishingly small for extremely short times, so that the effects caused by its vanishing 
are sufficiently negligible, as we have showed in Section |5.2| when we have analysed the 
gravitational potential and the UDM density contrast 5 = 5p/p&, Eq. (5.6). 

Studying observational constraints on our UDM models from SN 1A, CMB anisotropics 
and from the formation of the large-scale structure in the Universe will be the subject of 
a future analysis. In particular, it will be interesting to see if the feature wde < — 1 our 
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models show after the transition will be compatible with observations. On the theoretical 
side it will be important to check if, in a Universe filled with the scalar field component 



defined by the Lagrangian (A.l), the dynamics will display the desired behaviour, with the 
scalar field mimicking both the DM and DE components. Another issue concerning the 
dynamics of our UDM is the possible development of caustics in the non-linear regime. For 
instance, it is well known that inhomogeneous tachyon matter fluctuations could develop 
caustics, e.g. see |78], [79|. However, comparing with these models our Lagrangian has an 
additive potential V(4>). These open theoretical issues will be analysed in a forthcoming 
work. 
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A. Explicit reconstruction of Lagrangian C(<p, Y) 

Following the prescriptions in Section |3|, the analytic expression of the scalar field La- 
grangian for our UDM model turns out to be 



£«>, Y) = -f(4>) 



1 - 2 



K4>) 



i + 



3 <j> 2 



Y - V(d>) 



(A.l) 



Now, using the parametric equation of state p(N)-p(N) defined by Eqs. ([L4|)- (|4.5| ) and the 
speed of sound defined in Eq. ( |5.1[) , we can infer a relation among the parameters /i, v, Cqo 
and /3, i.e. 



n = Mi - 4>) 



(A.2) 



- 21 - 



Thanks to this relation we can eliminate fi, finally obtaining the following potentials 



1 ~ 1 

— — tanh 



1/2 



1 ' I 
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where cp t = [2/ (3z^a t " 3 )] 
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